Claims 




1. ^\ method of docking a ligand to a protein 
comprising: 

per^forming a pre-docking conformational search to 
generate Viultiple solution conformations of the ligand; 
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generating a binding site image of the protein, 
said bindin\g site image comprising multiple hot spots; 
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matchinlg hot spots of the binding site image to 
atoms in at least one solution conformation of the 
multiple solution conformations of the ligand to obtain 
at least one Ligand position relative to the protein in 
a ligand-protiin complex formation; and 

optimizing the at least one ligand position while 
allowing translation, orientation and rotatable bonds 
of the ligand td vary, and while holding the protein 
fixed. > 

2. The metW)d of claim 1, wherein said performing the 
pre-docking conf oriiational search comprises creating a 
database of the multiple solution conformations and storing 
said three-dimensionaV database for subsequent use by said 
matching. 



1 
2 
3 



3. The method of c^laim 2, wherein said database of 
the multiple solution con^rmations comprises a 
conformational database of V combinatorial library. 
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4. TUie method of claim 1, wherein said performing the 
pre-docking\conformational search comprises: 

randomly generating a plurality of uniformly 

4 distributled conformations of the ligand; 

5 minimizing a strain of each conformation of the 

6 plurality oVf uniformly distributed conformations; 

7 using t|he strain and a solvent accessible surface 

8 area of each^conformation to rank the conformations; 

9 and 

I 

'fl clusterirb the conformations and retaining a 

:i| desired number\of top clusters of conformations, said 

B ^2 retained number! of top clusters of conformations 

comprising said\multiple solution conformations of the 



14 ligand. 



1 5. The method of claim 1, wherein said generating the 

binding site image includes at least one of creating a list 

3 of apolar hot spots identifying points in the binding site 

4 that are favorable for an apolar atom to bind, and 

5 generating a list of polar hot spots identifying points in 

6 the binding site that are favorable for a hydrogen bond 

7 donor or acceptor to bind. 



6. The mVthod of claim 5, wherein said generating the 
binding site imaye further comprises: 

placing^ grid around the binding site of the 
protein; 
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grid; 



deteWining a hot spot search volume using said 



determining hot spots using a grid-like search of 
the hot spotXsearch volume; and 
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for each Vype of hot spot, clustering the hot 
spots and retailing a desired number of clusters of hot 
spots with bestWores, said desired number of clusters 
comprising said i)i^ultiple hot spots to be employed by 
said matching. 
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7. Th^ method of claim 1, wherein said matching 
comprises : 



matching atoms of the at least one solution 
'conformation to appropriate hot spots of the protein by 
positioning Vhe at least one solution conformation as a 
rigid body in\o the binding site image; 

defining aXmatch, said match determining a unique 
rigid body transformation; and 
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using the uniVue rigid body transformation to 
place the at least Vne solution conformation of the 
ligand into the binAng site of the prot( 



;ein, 



. \ The method of claim 7, wherein sal 



d determining 



the uniqu^ rigid body transformation comprises determining 
(^3/ the unique\ rigid body transformation that minimi 



iimizes : 
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I(R, T) = |Hj - RAj - T 

j=i 



where 
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E, = 



R = 
T = 



a j^"^ ftot spot of the protein; 



a atom of the at least one solution 
con forme t ion; 

a 3x3 rotation matrix; and 
a translation vector. 
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9. The method oflclaim 1, wherein said optimizing 
comprises optimizing mulltiple protein-ligand complex 
formations, said optimizing comprising: 
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eliminating each ligand position having a 
predetermined percenjiage of ligand atoms with a steric 
clash; 
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ranking remaining ligand positions using an atom 
pairwise score with a Idesired atom score cutoff; 

after ranking, clustering the ligand positions and 
selecting a top number \ of ligand positions; and 

optimizing each ligabd position of the n 
positions, allowing the translation, rotation and 
rotatable bonds of the lig\nd to vary. 
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10. The method of claim 9, wherein said optimizing 
comprises optimising each ligand position of the n positions 
using a BFGS optiUzation algorithm with a simple atom 
pairwise score, alMwing the translation, rotation and 
rotatable bonds of tVe ligand to vary. 
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11. A system for docking a ligand to a prot 
comprising: 



em 
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^ans fbr performing a pre-docking conformational 
search to gei^rate multiple solution conformations of 
the ligand; 



6 
7 



means for 
protein, said 
hot spots; 



generating a binding site image of the 
'inding site image comprising multiple 



rt 
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means for matching hot spots of the binding site 
image to atoms kn at least one solution conformation of 
the multiple solution conformations of the ligand to 
obtain at least |one ligand position relative to the 
protein; and 
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means for ojltimizing the at least one ligand 
position while alkowing translation, orientation and 
rotatable bonds ot the ligand to vary, and while 
holding the protean fixed. 

12. Tfte system of claim 11, wherein said means for 
performing thV pre-docking conformational search comprises 
means for creating a database of the multiple solution 
conformations aXd for storing said three-dimensional 
database for sub^^uent use by said matching. 

13. The system of claim 12, wherein said database of 
the multiple solutioA conformations comprises a 
conformational databak of a combinatorial library. 
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system of claim 11, wherein said means for 
performing th^ pre-docking conformational search comprises: 

meanfc for randomly generating a plurality of 
uniformly distributed conformations of the ligand; 

means Ifor minimizing a strain of each conformation 
of the plurklity of uniformly distributed 
conformations; 
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means rior using the strain and a solvent 
accessible sUface area of each conformation to rank 
the conformations; and 

means foW clustering the conformations and 
retaining a dLired number of top clusters of 
conformationsA said retained number of top clusters of 
conformations domprising said multiple solution 
conformations o\ the ligand. 

15. The system of claim 11, wherein said means for 
generating the binding site image includes at least one of 
means for creating a list of apolar hot spots identifying 
points in the binding site that are favorable for an apolar 
atom to bind, and means for generating a list of polar hot 
spots identifying points in the binding site that are 
favorable for a hydrogen bond donor or acceptor to bind. 



r \b\ ' V' °^ ^^^^"^ ^h^^^i^ ^^id means for 

9^1^ generating Vhe binding site image further comprises: 
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means Ifor placing a grid ar 
of the protein; 



ound the binding site 
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means f^r determining a hot spot search vol 
using said griid; 
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means fori determining hot spots using a grid-like 
search of the lot spot search volume; and 

for each tjrpe of hot spot, means for clustering 
the hot spots ark for retaining a desired number of 
clusters of hot ^.ots with best scores, said desired 
number of cluster\^ comprising said multiple hot spots 
to be employed by ^aid matching. 

17. \The system of claim 11, wherein said means for 
.matching opmprises: 



m^ans for matching atoms of the at least one 
solutioA conformation to appropriate hot spots of the 
protein W positioning the at least one solution 
conformat\^on as a rigid body into the binding site 
image; 
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means flor defining a match, said match determining 
a unique rigi\l body transformation; and 

means for Vising the unique rigid body 
transformation Ao place the at least one solution 
conformation of \he ligand into the bindi 



protein. 



ng site of the 
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18. The ^ystem of claim 17, wherein said determining 
the unxque rigik body transformation comprises determining 
the unique rigici^ body transformation that minimizes- 
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where : 



R = 
T = 



3 

IKR, T) = J |Hj - RAj - tI 

j=i 



hot spot of the protein; 



a j^^^tom of the at least one solution 



conformation; 
a 3x3 dotation matrix; and 
a translation vector. 
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19. The system o\ claim 11, wherein said means for 
optimizing comprises m^ans for optimizing multiple protein- 



ligand complex formati 
comprising: 



ohs, said means for optimizing 
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means for eliminating each ligand position having 
a predetermined percentage of ligand atoms with a 
steric clash; 
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means for rankirig remaining ligand positions using 
an atom pairwise scor^ with a desired atom score 
cutoff; 
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afteV ranking, means for clustering the ligand 
positions Wd selecting a top number n of ligand 
positions; Vnd 



Q 
iJl 
111 



14 
15 
16 

1 
2 
3 
4 
5 
6 



means foV optimizing each ligand position of the n 
positions, allUing the translation, rotation and 
rotatable bond^ of the ligand to vary. 

20. The systemW claim 19, wherein said means for 
optimizing comprises leans for optimizing each ligand 
position of the n posiVions using a BFGS optimization 
algorithm with a simplA atom pairwise score, allowing the 
translation, rotation aAd rotatable bonds of the ligand to 
vary. 
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21. At least one program storage device readable by a 
machine, tangiky embodying at least one program of 
instructions exkutable by the machine to perform a method 
of docking a ligVnd to a protein, comprising: 



perfo 
generate mul 



rir^ing a pre-docking conformational search to 
iple solution conformations of the ligand; 



generatjng a binding site image of the protein, 

site image comprising multiple hot spots; 



said binding 

matching 
atoms in at 1 
multiple solu 

at least one |.igand position relative to the protein; 
and 



hot spots of the binding site image to 
east one solution conformation of the 
tion conformations of the ligand to obtain 



mg 



optimizi 
allowing trans 
of the ligand 
fixed. 



the at least one ligand position while 
ation, orientation and rotatable bonds 
\o vary, and while holding the protein 
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22. ke at least one program storage device of claim 
21, wherein^aid performing the pre-docking conformational 
search comprises creating a database of the multiple 
solution conformations and storing said three-dimensional 
database for subsequent use by said matching. 
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23. The at. least one program storage device of cl 
22, wherein said\database of the multiple solution 
conformations comprises a conformational database of a 
combinatorial library. 
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at least one program storage device of claim 
21, wherein\said performing the pre-docking conformational 
search comprVses: 



ranjomly generating a plurality of uniformly 
distribuljed conformations of the ligand; 

miniiiizing a strain and a solvent accessible 
surface aiea of each conformation of the plurality of 
uniformly jflistributed conformations; 
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using the strain of each conformation to rank the 
conformatioiis/ and 

clustering the conformations and retaining a 
desired numbir of top clusters of conformations, said 
retained numller of top clusters of conformations 
comprising saljd multiple solution conformations of the 
ligand. 

25. The at least one program storage device of claim 
21, wherein said generating the binding site image includes 
at least one of creating a list of apolar hot spots 
Identifying points in the binding site that are favorable 
for an apolar atom to bind, and generating a list of polar 
hot spots identifying points in the binding site that are 
favorable for a hydrogen bond donor or acceptor to bind 



1 26. \ The at least one program storage device of claim 

2 25, whereiyn said generating the binding site image further 
comprises 
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placing a grid around the binding site of the 
protein; 
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deteiWining a hot spot search volume using said 



grid; 



deterriining hot spots using a grid-like search of 
the hot spot search volume; and 
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for eack type of hot spot, clustering the hot 
spots and retWning a desired number of clusters of hot 
spots with belt scores, said desired number of clusters 
comprising saijl multiple hot spots to be employed by 
said matching. 



27. The\at least one program storage device of cla 
21, wherein sa\d matching comprises: 



im 



matching atoms of the at least one solution 
conformation! to appropriate hot spots of the protein by 
positioning he at least one solution conformation as a 
rigid body in\o the binding site image; 

defining aUatch, said match determining a unique 
rigid body transformation; and 

using the unLue rigid body transformation to 
place the at least\one solution conformation of the 
ligand into the binding site of the protein. 
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ThV at least one program storage device of claim 



^^2/> 27, wherein sWd determining the unique rigid body 
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transformati 



transformatilpn comprises determining the unique rigid body 
)n that minimizes: 



where : 



I(R, T) = 2] |h. - RAj - t| 
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T = 



hot spot of the protein; 



a t-th 

a ^tom of the at least one soluti 
conVformation; 

a 3^3 rotation matrix; and 
a trl^nslation vector. 
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29. The at lea^t one program storage device of claim 
21, wherein said optimizing comprises optimizing multiple 
protem-ligand comple)^ formations, said optimizing 
comprising: 
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eliminating ekch ligand position having a 

predetermined perce^age of ligand atoms with a steric 
clash; 
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11 



ranking remaining ligand positions using an atom 
pairwise score with a ijesired atom score cutoff; 



after ranking, clu£ 
selecting a top number n 



tering the ligand positions and 
of ligand positions; and 
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optimizing 



each ligand position of the n 
positions, alloWng the translation, rotation and 
rotatable bonds Vf the ligand to vary. 

30. The at leastWe program storage device of claim 
29, Wherein said optimikng comprises optimizing each ligand 
position of the n positions using a BFGS optimization 
algorithm with a simple itom pairwise score, allowing the 
translation, rotation and\rotatable bonds of the ligand to 
vary. 
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